
EFFECTS OF ENDWALL GEOMETRY AND STACKING ON 
TWO-STAGE SUPERSONIC TURBINE PERFORMANCE 


Daniel J. Dorney* Lisa W. Griffin 1 

Fluids Dynamics Analysis Branch 
NASA Marshall Space Flight Center 
Marshall Space Flight Center, Alabama, USA 

Frank W. Huber* 

Riverbend Design Services 
Palm Beach Gardens, Florida, USA 

Douglas L. Sondak§ 

Office of Information Technology 
Boston University 
Boston, MA, USA 


ABSTRACT 


NOMENCLATURE 


The drive towards high-work turbines has led to 
designs which can be compact, transonic, super- 
sonic, counter rotating, or use a dense drive gas. 
These aggressive designs can lead to strong sec- 
ondary flows and airfoil flow separation. In many 
cases the secondary and separated flows can be 
minimized by contouring the hub/shroud endwalls 
and/or modifying the airfoil stacking. In this study, 
three-dimensional unsteady Navier-Stokes simula- 
tions were performed to study three different end- 
wall shapes between the first-stage vanes and rotors, 
as well as two different stackings for the first-stage 
vanes. The predicted results indicate that chang- 
ing the stacking of the first-stage vanes can sig- 
nificantly impact endwall separation (and turbine 
performance) in regions where the endwall profile 
changes. 


C - Axial chord 

/ - Frequency 

M - Mach number 

P - Static pressure 

r - Radius 

A P - Amplitude of pressure variation 

f - (r - V hub) /{ r tip - Thub) 

T - Static temperature 

W - Work 


GREEK SYMBOLS 

fj - Relative circumferential angle 
7] - Efficiency 

7 - ratio of specific heats 

p - Density 

0 - Rotor rotational speed 
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SUBSCRIPTS 

in - Inlet 

out - Outlet 

t - Stagnation quantity, time derivative 
ts - Total-to-static 

tt - Total-to-total 
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0 - Vane-1 inlet 

4 - Rotor-2 exit 

1/2 - Half-amplitude 

oo - Free stream 

INTRODUCTION 

Modern high- work turbines can be compact, tran- 
sonic, supersonic, counter rotating, or uses a dense 
drive gas. The vast majority of modern rocket tur- 
bine designs fall into these categories. These tur- 
bines are often characterized by large amounts of 
flow unsteadiness. The flow unsteadiness can have 
a major impact on the turbine performance and 
durability. For example, the Space Transportation 
Main Engine (STME) fuel turbine, a high- work tran- 
sonic design, was found to have an unsteady interrow 
shock which reduced efficiency by 2 points and in- 
creased dynamic loading by 24 percent. The Revolu- 
tionary Reusable Technology Turbopump (RRTT), 
which uses full flow oxygen for its drive gas, was 
found to shed vortices with such energy as to raise 
serious blade durability concerns. In both cases, the 
sources of the problems were uncovered (before tur- 
bopump testing) with the application of validated, 
unsteady computational fluid dynamics (CFD) to 
the designs. In the case of the RRTT and the Al- 
ternate Turbopump Development (ATD) turbines, 
the unsteady CFD codes were used not just to iden- 
tify problems, but to guide designs which mitigate 
problems due to unsteadiness. Using unsteady flow 
analyses as a part of the design process has led to 
turbine designs with higher performance (which ef- 
fects temperature and mass flow rate) and fewer dy- 
namics problems. The works of Griffin et al. [l]-[4], 
Garcia et al. [5] and Griffin and Dorney [6] are ex- 
amples of the application of unsteady CFD to rocket 
turbine designs. 

More recently, CFD has been used to design a 
two-stage supersonic turbine which will be tested 
experimentally during 2002 [7]-[9]. Numerical sim- 
ulations (including meanline, two-dimensional CFD 
and three-dimensional CFD analyses in conjunction 
with optimization techniques) were used to design 
both the flowpath and the airfoil geometries. Dur- 
ing the course of this work a large separated flow 
region was detected on the hub endwall between the 
first-stage vane and the first-stage rotor. 

Two methods used are normally to control the 
secondary / separated flows (and associated losses) in 
supersonic turbines: endwall contouring and airfoil 
stacking. In the current investigation the flow path 
between the first-stage vanes and rotors, and the 


stacking of the first-stage vanes, were varied in an ef- 
fort to improve turbine performance. The geometric 
variations have been studied by performing a series 
of unsteady three-dimensional numerical simulations 
for the two-stage turbine. 

NUMERICAL ALGORITHM 

The governing equations considered in this 
study are the time dependent, three-dimensional 
Reynolds-averaged Navier-Stokes equations. To ex- 
tend the equations of motion to turbulent flows, an 
eddy viscosity formulation is used. The turbulent 
viscosity is calculated using the two-layer Baldwin- 
Lomax algebraic turbulence model [10]. 

The numerical algorithm used in the three- 
dimensional computational procedure consists of 
a time-marching, implicit, finite-difference scheme. 
The procedure is third-order spatially accurate and 
second-order temporally accurate. The inviscid 
fluxes are discretized according to the scheme devel- 
oped by Roe [11]. The viscous fluxes are calculated 
using standard central differences. An approximate- 
factorization technique is used to compute the time 
rate changes in the primary variables. Newton sub- 
iterations are used at each global time step to in- 
crease stability and reduce linearization errors. For 
all cases investigated in this study, two Newton sub- 
iterations were performed at each time step. Mes- 
sage Passing Interface (MPI) and OpenMP software 
have been implemented into the numerical analysis 
to reduce the computation time for large-scale three- 
dimensional simulations. 

The Navier-Stokes analysis uses O- and H-type 
zonal grids to discretize the flow field and facili- 
tate relative motion of the rotating components (see 
Fig. 1). The O-grids are body-fitted to the surfaces 
of the airfoils and generated using an elliptic equa- 
tion solution procedure. They are used to prop- 
erly resolve the viscous flow in the blade passages 
and to easily apply the algebraic turbulence model. 
The algebraically-generated H-grids are used to dis- 
cretize the remainder of the flow field. 

The computational analysis has been validated on 
several supersonic turbine geometries ( e.g ., Refs. [6], 
[12], [13]). 

BOUNDARY CONDITIONS 

The theory of characteristics is used to determine 
the boundary conditions at the inlet and exit of the 
computational domain. For subsonic inlet flow four 
quantities are specified and one is extrapolated from 
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the interior of the computational domain. In partic- 
ular, the total pressure, total temperature, and the 
circumferential and radial flow angles are specified 
as a function of the radius. The upstream running 
Riemann invariant is extrapolated from the interior 
of the computational domain. 

For subsonic outflow one flow quantity is speci- 
fied and four are extrapolated from the interior of 
the computational domain. The circumferential and 
radial flow angles, total pressure, and the total tem- 
perature are extrapolated from the interior of the 
computational domain. The pressure ratio, Pa/Ph* 
is specified at mid-span of the computational exit 
and the pressure at all other radial locations at the 
exit is obtained by integrating the equation for ra- 
dial equilibrium. For supersonic outflow all the flow 
variables are extrapolated. Periodicity is enforced 
along the outer boundaries of the H-grids in the cir- 
cumferential direction. 

For viscous simulations, no-slip boundary condi- 
tions are enforced along the solid surfaces. It is as- 
sumed that the normal derivative of the pressure is 
zero at solid wall surfaces. In addition, a specified 
heat flux distribution is held constant in time along 
the solid surfaces. 

The flow variables at zonal boundaries are explic- 
itly updated after each time step by interpolating 
values from the adjacent grid. 

GEOMETRY AND GRIDS 

The two-stage supersonic turbine configuration, 
typical of those proposed for a reusable launch ve- 
hicle, has 12 first-stage vanes, 30 first-stage ro- 
tors, 73 second-stage vanes and 56 second-stage ro- 
tors. In the current effort, a 1 5- vane/ 30- rotor / 75- 
vane/60-rotor (1/2/5/4) airfoil approximation has 
been made. To keep the pitch-to-chord ratio (block- 
age) constant, the first-stage vanes were scaled by a 
factor of 12/15, the second-stage vanes were scaled 
by a factor of 73/75 and the second-stage rotors were 
scaled by a factor of 56/60. The tip clearance in the 
first- and second-stage rotors was set at the design 
value of approximately 2.0% of the respective rotor 
heights. 

The grid densities (number of passagesxixjxk) 
for the turbine simulations are presented in Table 1. 
The total number of grid points used to discretize 
the turbine was 4,139,957. Figure 1 illustrates an 
(x - y) view of the grids at midspan, where every 
other grid point in each coordinate direction has 
been removed for clarity. Figure 2 illustrates the 
grids used to discretize the clearance region of the 


second-stage rotor. Figure 3 shows a perspective 
view of the two-stage turbine. The average value 
of , the non-dimensional distance of the first grid 
line above the surface was approximately 1.0 for the 
airfoils surfaces and 1.5 for the endwall surfaces. 

The simulations were run on 24 (400 MHz) pro- 
cessors of an SGI Origin2000. Each simulation was 
run for 15.0 global cycles (one complete rotor revo- 
lution) at 22,000 iterations per cycle. A global cycle 
is defined as the time it takes for the two first-stage 
rotor blades to pass by the first-stage vane airfoil. 
The value of 22,000 iterations per cycle was chosen 
to resolve all the (expected) frequencies of interest. 
Each iteration required approximately 9.0 seconds 
computation time on 24 processors. The time peri- 
odicity of the solutions was determined by interro- 
gating pressure traces at different points along the 
airfoil surfaces. 

NUMERICAL RESULTS 

The two-stage turbine under consideration has a 
design inlet Mach number of Mo = 0.08, an inlet 
static pressure of 2225 psia , and an inlet static tem- 
perature of approximately To = 2225 R. The rotor 
rotates at Q = 31,343 RPM, the Reynolds number 
(based on the inlet conditions and the rotor axial 
chord) is approximately 1.2 x 10 6 and the ratio of the 
rotor exit static pressure to vane inlet total pressure 
is Pa/Pio = 0.1135. The operating fluid is hydrogen- 
rich steam and the average ratio of specific heats is 
7 = 1.3538. 

Four different simulations have been performed to 
determine the effects of endwall shape and first-stage 
vane stacking on the performance of the two-stage 
turbine. 

1. Case 1 - The radius of the inner diameter (ID) 
and outer diameter (OD) are transitioned from 
the nominal values in the first-stage vane pas- 
sage to the final values in the rotor passage be- 
ginning three-quarters of the way between the 
vane trailing edge and concluding at the rotor 
leading edge (see Fig. 4). Case 1 represents 
the original turbine geometry based on previ- 
ous work [7]- [9]. 

2. Case 2 - The radius of the ID and OD are tran- 
sitioned from the nominal values in the vane 
passage to the final values in the rotor passage 
beginning at the vane trailing edge and conclud- 
ing at the rotor leading edge (see Fig. 4). 

3. Case 3 - The radius of the inner diameter (ID) 
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is kept constant at the value used in the rotor. 
The radius of the OD is transitioned from the 
nominal value in the vane passage to the final 
value in the rotor passage beginning at the vane 
trailing edge and concluding at the rotor leading 
edge (see Fig. 4). Note, in this case the height 
of the vane was increased to keep the vane flow 
area the same as in Cases 1 and 2. 

4. Case 4 - The flow path is identical to that used 
in Case 3. The vane airfoils are stacked along 
the trailing edge instead of the center of gravity 
(as it was in Cases 1, 2 and 3). 

An overview of the turbine flow field is presented in 
Figs. 5 and 6, which contain instantaneous absolute 
Mach and entropy contours, respectively, at midspan 
of the turbine for Case 4. Figure 5 highlights the 
strong expansion and shock structures between the 
first-stage vane and rotors, as well as weaker shock 
structures near the leading and trailing edges of the 
second-stage vane. Figure 6 highlights flow separa- 
tion on the suction surface of both rotors, which is 
typical of supersonic turbine flow fields. 

Time-averaged entropy contours midway between 
the first-stage vane and rotor passages for Case 1 
are shown in Fig. 7. The contours show a large re- 
gion of separated flow extending from the hub to 
approximately 20% of the span. The axial extent 
of the region was confined to the area between the 
vane trailing edge and the rotor leading edge. It was 
initially theorized that the large separated flow re- 
gion was being induced by the rapid expansion in 
the endwall flowpath in Case 1. Reducing the slope 
of the endwall in Cases 2 and 3, however, did not sig- 
nificantly affect the size of the separated flow region 
(see Figs. 8 and 9). The next hypothesis for the large 
separated flow region was that stacking the vanes 
along the center of gravity causes the vane throat to 
point outwards towards the shroud endwall, giving 
the flow a tendency to pull away from the hub end- 
wall. Re-stacking the vanes along a radial line con- 
necting the trailing edge points significantly reduced 
the size of the separated flow region (see Fig. 10). 

Tables 2 to 5 contain the time-averaged relative- 
frame flow quantities at the inlet and exit of each 
blade row for all four cases. Some of the relevant 
information which can be deduced from these tables 
includes: 

• Reducing the size of the separated flow region 
in Case 4 resulted in a value of the average vane 
exit Mach number (1.37) which is closer to the 
design value of 1.50. 


• Reducing the size of the separated flow region 
in Case 4 gives a significant increase (nearly 6 
points) in turbine efficiency compared to Case 
1. A more detailed comparison of Cases 1 and 
4 is presented below. 

• The changes made to the first stage had little 
effect on the flow in the second stage. 

The differences between the flow fields in Cases 1 
and 4 are explored in more detail by comparing sur- 
face pressures and radial profiles. Figures 11 and 12 
illustrate the time-averaged surface pressures on the 
first-stage vane and rotor, respectively, at 10%, 50% 
and 90% of the span. Re-stacking the vane results 
in the loading shifting towards the leading edge at 
the hub and towards the trailing edge at the tip (see 
Fig. 11). The vane loadings are similar at 50% of the 
span. There are significant differences between the 
rotor surface pressures in Cases 1 and 4 (see Fig. 12). 
At 10% span the surface pressures are much greater 
in Case 1 because of the low Mach number associ- 
ated with the separated flow. The loading on the 
rotor at 10% span is also greater in the Case 1 than 
in Case 4. The presence of the low flow region near 
the hub in Case 1 forces more of the flow through 
the outboard regions of the rotor. The loadings at 
50% and 90% span are similar in Cases 1 and 4. 

An indication of the unsteadiness in the turbine 
can be obtained from Figs. 13-16, which show the 
unsteady pressure envelopes for all four blade rows 
at 10%, 50% and 90% of the span in Case 4. The un- 
steadiness on the first-stage vane, which is generated 
mainly by interaction with the rotor potential field 
and bow shock, is confined to the unconvered por- 
tion of the suction surface by the choked throat in 
passage (see Fig. 13). The first-stage rotor blades, as 
well as the airfoils in the second stage of the turbine, 
experience a significant amount of unsteadiness. The 
unsteadiness in these blade rows is generated from 
several sources, including interaction with the po- 
tential fields of the upstream and downstream blade 
rows, periodic interaction with the wakes from up- 
stream blade rows and temporal variations in shock 
location and strength. The unsteadiness in all four 
blade rows is relatively constant as a function of the 
span. 

Fourier decompositions of the unsteady pressure 
at 5% of the axial chord on the suction surface of 
the first-stage rotor at 10% and 50% of the span are 
shown in Figs. 17 and 18, respectively, for Cases 1 
and 4. Note, the frequencies in Case 1 were shifted 
by 1 kHz to facilitate comparisons. At 10% span the 
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reduction in the separated flow region causes the ro- 
tor to experience more influence from the vane wake 
(as indicated by an increase in the unsteadiness near 
8 kHz). At 50% span both cases exhibit significant 
unsteadiness at the vane passing frequency and twice 
the vane passing frequency. The higher harmonic 
may actually be caused by the vane trailing edge ex- 
pansion fan, which is offset from the vane wake by 
nearly half the vane circumferential pitch. 

Time-averaged radial profiles of the absolute Mach 
number at the exit of each blade row in Cases 1 and 
4 are shown in Figs. 19-22. The Mach number profile 
at the exit of the first-stage vane (see Fig. 19) clearly 
shows the extent of the separated flow region. The 
flow field begins to recover by the time it exits the 
first-stage rotor (see Fig. 20), and the Mach number 
profiles are nearly identical in the second stage of 
the turbine (see Figs. 21 and 22). The presence of 
tip leakage flow is evident behind both rotor rows 
(see Figs. 20 and 22). 

Time-averaged radial profiles of the absolute cir- 
cumferential flow angle at the exit of each blade row 
in Cases 1 and 4 are shown in Figs. 23-26. As noted 
above, the large flow (and angle) deficit near the 
first-stage vane hub in Case 1 (see Fig. 23) causes a 
radial shift in the flow distribution within the rotor 
(see Fig. 24). Similar to the Mach number profiles, 
the flow angle profiles in the second stage are sim- 
ilar for both cases. Thus, the geomertic variations 
in the first-stage do not have a significant impact on 
the flow in the second stage. 

For completeness, the time- averaged radial profiles 
of the absolute total pressure are shown in Figs. 27- 
30. As expected, the character of the total pressure 
profiles follows closely with that of the Mach number 
profiles. The large separated flow region at the exit 
of the first-stage vane, the redistribution of the flow 
in the first-stage rotor, as well as the recovery of the 
flow in the second stage are all visible. 

CONCLUSIONS 

A set of unsteady three-dimensional Navier-Stokes 
simulations has been used to investigate the effects 
of endwall shape and first-stage vane stacking on the 
performance of a two-stage supersonic turbine. Re- 
stacking of the vanes was successfully used to elimi- 
nate a large separated/secondary flow region at the 
hub between the first-stage vanes and rotors. Alter- 
ing the shape of the endwall in the first stage had 
little effect on the separated flow region. There was a 
significant performance increase obtained at the de- 
sign flow conditions by reducing the separated flow 


region. It is anticipated that the benefits of improv- 
ing the behavior of the flow near the endwall will be 
even greater at off-design operating conditions. 
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Grid Type 


Vane-1 


Rotor- 1 


Vane-2 


Rotor-2 


Tip 

Total Points 


1x141x31x51 

1x67x41x51 


363,018 


2x151x21x51 

2x64x41x51 

2x151x26x7 

646,054 


5x121x31x51 

5x64x41x51 

1,525,625 


4x121x31x51 

4x82x41x51 

4x121x16x7 

1,505,260 


Table 1: Grid dimensions for the 2-stage turbine. 


Variable 

Vane-1 

Rotor- 1 

Vane-2 

Rotor-2 

M in 

0.08 

0.91 

0.59 

1.09 

Mout 

1.17 

0.85 

1.09 

0.81 

P in (psia) 

2225. 

604. 

564. 


Pout (psia) 

mm i 

564. 



Pt in (psia) 

2235. 



458. 

Rt otl f (psia) 

1611. 

923. 



Tt in (R) 

2225. 

2000. 

1863. 


Tt out (R) 

2224. 

1986. 


PVfrWIflM 

0 in (deg) 

0.0 

60.5 

-57.1 

56.8 

Pout (deg) 

63.8 

-70.1 

67.4 

-55.2 

W (BTU/lbm) 


495. 

- 

418. 

Reaction 

- 


- 


T)tt — overall 

- 

- 

- 


Tits— overall 

- 

- 

- 

wmmm 


Table 2: Turbine time-averaged Sow quantities for Case 1. 


Variable 

M in 

Mout 

Pin (psia) 

Pout (psia) 

PUn (psia) 
Ptput (psia) 
Ttjn (R) 

Ttout (R) 

Pin (deg) 

Pout (deg) 

W (BTU/lbm) 
Reaction 


Vane-1 

008 

1.15 

2225. 

611. 

2235. 

1593. 

2225. 

2224. 

0.0 

61.0 


tyt — o verall 
rjta— overall 


Rotor- 1 

090 

0.86 

611. 

560. 

1124. 

922. 

2004. 

1989. 

59.2 

-70.7 

502. 

0.066 


Vane-2 

059 

1.09 
560. 
299. 
719. 
626. 
1865. 
1856. 
-56.5 

66.9 


Rotor- 2 
T09 
0.80 
299. 

254. 

455. 

408. 

1706. 

1705. 

56.8 

-55.1 

416. 

0.159 

0.623 

0.568 


Table 3: Turbine time-averaged Sow quantities for Case 2. 
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Variable 

Vane-1 

Rotor- 1 

Vane-2 

Rotor-2 

M in 

0.08 


0.59 

0.79 

Mout 

1.15 

0.86 

1.08 

0.85 

P in (psia) 

2225. 

598. 



Pout (psia) 

598. 

554. 



Ptin (psia) 

2235. 

1115. 

711. 

m 

P^out (psia) 

1577. 

922. 

622. 

406. 

Tt in (R) 

2225. 




Ttout (R) 

2224. 

1989. 

1857. 


0in (deg) 

0.0 

59.1 

-57.3 

56.7 

Pout (deg) 

60.0 

-70.7 

67.5 

-55.3 

W (BTU/lbm) 

- 

481. 

- 

425. 

Reaction 

- 

0.059 

- 

0.165 

'Hit — overall 

- 

- 

_ 

0.609 

qt8— overall 

- 

- 

- 

0.553 


Table 4: Turbine time-averaged Sow quantities for Case 3. 


Variable 

Vane-1 

Rotor- 1 

Vane-2 

Rotor-2 

Mi n 

0.08 


0.59 

0.80 

M 0 ut 

1.37 

0.86 

1.09 

0.86 

P in (psia) 





Pout (psia) 

602. 

560. 

300. 


Ptin (psia) 

2235. 

1237. 


454. 

Ptout (psia) 

1835. 

926. 

625. 

408. 

Tt in (R) 


1996. 

1863. 

1707. 

Ttout (R) 

2223. 

1989. 

1856. 

1703. 

a in (deg) 

0.0 

72.7 


56.8 

&out (deg) 

76.4 

-69.7 


-55.1 

W (BTU/lbm) 

- 

583. 

- 

427. 

Reaction 

- 


- 

0.164 

'Htt — overall 

- 


- 

0.680 

'Hta — overall 

- 

- 

- 

0.618 


Table 5: Turbine time-averaged Sow quantities for Case 4. 
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Figure 10: Time-averaged entropy contours - vane-1 
exit - Case 4. 
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Figure 11: Time-averaged pressure - vane-1. 
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Figure 12: Time-averaged pressure - rotor-1. 


Figure 13: Pressure envelope - Case 4 - vane-1 
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Figure 14: Pressure envelope - Case 4 - rotor-1. Figure 15: Pressure envelope - Case 4 - vane-2. 
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Figure 16: Pressure envelope - Case 4 - rotor-2. 


Figure 18: Fourier decomposition of unsteady pres- 
sure - rotor-1 - 5% axial chord - suction surface - 
50% span. 


14 













Figure 24: Absolute circumferential flow angle pro- 
file - rotor- 1 exit. 



Figure 26: Absolute circumferential How angle pro- 
file - rotor-2 exit. 











